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MOLECULAR  THEORY  OF  SOLUTIONS  AND  EXTRACTION 


OF  METALS  AND  ACIDS  BY  ASSOCIATED  REAGENTS 


Ye.  V.  Komarov 


During  the  investigation  of  extraction  of  actinide  and  other 
elements  and  acids  by  such  associated  reagents  as  organophosphorous 
acids,  the  salts  of  amines  and  others,  frequently  dependences  of 
the  coefficients  of  distribution  of  substance  on  the  concentration 
of  the  extracting  reagent  are  observed  which  are  very  simple  in 
overall  nature.  There  is  surprise  in  the  fact  that,  for  example, 
the  extraction  of  nitric  acid  from  an  aqueous  solution  of  nitrates 
of  amines  can  be  described  within  the  framework  of  the  simplest 
stoichiometric  relationships  without  taking  polymerization  into 
account  [1],  without  taking  into  account  the  strong  deviation  of 
the  behavior  of  such  solutions  from  the  laws  of  behavior  of  ideal 
solutions.  Just  as  incomprehensible  (from  the  point  of  view  of 
customary  concepts)  are  the  linear  dependences  (in  logarithmic 
scale)  of  the  coefficients  of  distribution  of  different  metals  on 
the  concentration  of  reagent  in  the  case  of  extraction  by  amines 
[2-8],  monoalkyl-  and  monoarylphosphoric  acids  [9-13]. 

Attempts  at  taking  polymerization  into  account  should  unavoid¬ 
ably  lead  to  consideration  of  the  sum  of  reactions  of  the  form 


i  •  M  +  A,  ||(Ari 

where  M  -  acid,  salt  dr  ion  of  the  metal,  -  r-mer  of  the  rea¬ 
gent,  Kir  -  equilibrium  constants  of  the  corresponding  reactions. 

In  this  case  the  expression  for  the  coefficient  of  distribution  M 
should  include  the  following  sum 

2'W  |Arl~-A<r. 

r.<  iM(r 

where  a^  -  activity  of  M  in  the  aqueous  phase,  [A^]  -  concentration 

of  r-mer  in  the  organic  phase,  yM  -  activity  coefficient  of  a 

ir 

complex  of  the  composition  M^Ap,  Yr  -  coefficients  of  activity  of 
the  r-mer  of  reagent  A. 

In  order  to  describe  exactly  the  concentration  dependence  of 
the  coefficient  of  distribution  of  particles  M,  it  is  necessary  to 
know,  in  addition  to  the  variables  the  dependence  on  gross- 

concentration  of  the  reagent  of  the  content  of  different  associates 
[A  ]  and  the  variables  -----  for  the  different  complexes  and 

*  1  A!  j , 

associates.  The  selection  of  specific  conditions  for  carrying  out 
the  experiment  makes  it  possible  at  least  to  find  the  concentra¬ 
tions  of  [Ar]  by  independent  methods  of  physicochemical  analysis. 
The  experimental  determination  of  the  coefficients  of  activity  of 
individual  associates  and  complexes  is  just  as  difficult  (or  even 
impossible)  a  problem  as  the  determination  of  the  absolute  activi¬ 
ties  of  the  individual  ions  of  the  electrolyte.  Therefore  for  the 
calculation  of  the  variable  _1^-  we  have  no  other  route  except 
theoretical  calculation. 

Selection  of  the  molecular-rstatistical  model.  Molecular 
theoi  ies  of  associated  solutions  developed  rapidly  especially  in 
the  fifties.  The  successes  in  this  area  were  reflected  in  a  number 
o:’  monographs  and  reviews  [13—18 ] .  In  analyzing  the  dependences 
of  thermodynamic  functions  on  the  composition  of  associated  solu- 


2. 


tions  of  different  types,  Rowlinsdn  [i4]  came  to  the  conclusion 
that  the  sign  and  magnitude  of  thermodynamic  functions  of  mixing 
in  systems  with  strong  directed  bonds  between  the  molecules  are 
determined  mainly  by  the  change  in  the  number  of  such  bonds  during 
mixing  of  ah  associated  liquid  with  other  liquids.  This  shows 
that  out  of  the  different  variants  of  molecular  theories  of 
associated  solutions  that  one  should  be  selected,  in  which  in  the 
construction  of  the  statistical  sum  of  the  molecular  mixture  in 
an  evident  form  the  formation  of  associates  and  solvates  is  taken 
into  consideration  [13,  19-21] .  Since  such  a  route  takes  into 
account  the  contributions  from  the  change  in  internal  energy  and 
the  orderliness  of  molecules  in  the  system,  then  corrections  for 
the  effects  of  a  second  order  of  smallness,  conditioned  by  the 
change  in  the  molecular  environment  of  associates  and  solvates  with 
the  change  in  the  composition  of  the  system,  can  be  found  on  the 
basis  of  the  theory  of  strictly  regular  solutions  in  a  zero  approx¬ 
imation  [16]. 

The  special  features  of  the  selected  molecular-statistical 
model  of  a  solution  are  reflected  precisely  in  the  following  record¬ 
ing  of  the  statistical  sum  of  a  molecular  mixture 


t  (.Vp,  .vT) 


A'T)]KU 

!) 


ur 


(i) 


where  the  indices  D  and  T  are  the  generalized  symbols  of  the  sort 
of  all  possible  particles  (including  associates  and  solvates)  in 
the  system,  Np  -  number  of  particles  of  sort  D,  g(ND,  NT)  -  number 
of  configurations  of  the  mixture,  -  statistical  sum  based  on 
the  inner  states  of  a  particle  of  sort  D,  U(ND,  NT)  -  configuration 
energy  of  the  system,  depending  on  the  number  and  class  of  contacts 
of  all  possible  particles  in  the  system. 

The  variable  g  can  be  calculated  by  using  Guggenheim's  formula 
[16],  if  the  molecular  complexes  are  either  cyclic  or  noncyclic. 

In  the  case  of  particles  of  an  arbitrary  form  a  generalized  Guggen- 
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helm  formula  should  be  used  [22].  The  variable  U  can  be  calculated 
in  the  following  manner. 

A  molecule  of  a  given  sort  D  in  a  general  case  has  surface 
sectors  which  are  not  equivalent  to  each  other  in  an  energy  respect. 
Assume  e  and  v  are  generalized  symbols  of  the  class  of  contact 
sectors  of  the  molecule,  so  that  e  can  take  specific  values  of  a, 
b,  c,  etc.,  just  as  v  (the  second  symbol  is  introduced  for  indi¬ 
cating  the  difference  of  classes  of  contact  sectors)* 


#Here  the  symbols  developed  by  Barker  are  applicable  [23]. 


The  share  of  contact  sectors  of  molecule  D  in  class  £  is 

•ava  A  t  A  J  It  I — I  .  imK  /A  "1  t  i  ^  La1  a  a  a*  a]  ! — \  .  4*  T  /* 


designated  by  the  symbol  .  Evidently  5 

The  variable  Uzv  denotes  the  energy  of  the  contact,  formed  by  the 
contact  sectors  of  classes  E  and  V  of  particles  of  the  same  or 
different  sorts,  The  overall  number  of  all  possible  contacts  will 
be  equal  to 


0  2*0*0- 
D 


Here  zqp  -  number  of  nearest  neighbors  of  a  molecule  of  sort  D, 


z  -  coordination  number  of  the  lattice. 


The  share  of  contact  sectors  of  class  £  of  all  the  particles 
in  the  solution  will  be  equal  to 


Vt-2- 


It  is  easy  to  show  that  the  overall  number  of  contacts  of  type 
(£—  V)  will  be  equal  to  QV£V,  ,  and  their  energy  QV^V^'U^. 

The  dependence  of  the  configuration  energy  on  composition  of 
the  solution  will  be  expressed  as: 


N 


m 


(.3) 
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The  Gibbs  thermodynamic  potential  for  a  condensed  system  can  be 
found  from  (1): 


Gcz—kT 


|taf(flD.  iVT)  +  2 Nd  In  +5 -  -t)J. 


(4) 


&■  1 


Since  U(Ndj  nij)  is  determined  by  the  relationships  (2)  and  (3),  and 
ing(Np,  Nt)  by  formula  [22] 

+  |)(rb-?D)  *b  to  (2^)  i.  /; 


* 

&?. 
Ijg^  > 

fe- : 


where  -  number  of  methods  of  arrangement  of  the  molecule  on 
the  quasi-lattice  points  of  the  solution  in  the  case  of  the  same 
fixed  statistical  segment,  0~D  -  coefficient  of  symmetry,  r  -  number 
of  statistical  segments  of  the  molecule,  then  the  chemical  poten¬ 
tial  of  a  particle  of  sort  D  is  determined  by  the  equation  (cyclic 
structures  of  the  associates  and  solvates  are  absent) 


lip  —  kT  | —In  ~  '{'d  +  In  JVD  —  2  ?p  (n  ^  gp^pH-  rb  (t  —  Q  ^  rp^p|  4 

*  D  d  i 

+  f  9d  2  * (*?*',  +  V?K‘  ~ 


I,  V 


(5) 


Depending  on  the  method  of  expression  of  the  composition  of 
the  solution  and  the  selection  of  the  standard  state  of  the  com¬ 
ponent,  expression  (5)  will  acquire  a  different  form.  In  the 
majority  of  works  on  extraction  they  use  the  molar  concentrations 
as  the  method  for  expressing  the  composition  of  phases.  Therefore, 
in  this  case  there  is  sense  in  expressing  the  content  of  particles 
in  the  system  in  moles  per  liter,  and  as  the  standard  state  to 
select  a  solution  which  possesses  the  properties  of  an  infinitely 


5 


diluted  solution.  Then  the  coefficient  of  activity  of  a  particle 
of  sort  D  can  he  calculated  Using  the  following  formula: 


I..  ■; ,3  =-  -  la  y-~-  :  r„  (  -ST  -  l )  In  r 

'  rS  —  'll/ 1» 

to  ll 

>  (I  “  *)  • + 1  ]  2 TF [ ■?  ( 1 - 1  f )  +  («■.  -  rf )  -  (r, r,  -  r*K*)J, 


(6) 


where  C|j  -  molar  concentration  of  particles  of  sort  D,  and  the 
index  S  denotes  the  sort  of  particles  which  form  the  solvent. 

Numerical  appraisals  using  formula  (6)  are  possible  for 
specific  systems  and  specific  conditions  of  the  experiment.  Here 
we  will  consider  only  two  important  examples:  a  solution  of 
monooctylphosphoric  acid  in  CCl^  and  a  solution  of  trioctylammonia 
nitrate  in  benzene. 

Numerical  appraisals,  ^he  system  mono-N-octylphosphoric 
acid  (A)  -  complex  (B)  -  CO  .).  Since  an  ion  of  the  metal  in 
its  dimensions  is  considerably  smaller  than  the  dimensions  of 
molecules  of  A  and  S,  and  the  dimensions  of  A  are  approximately 
twice  the  dimensions  of  S,  then  rA  =rM^r=i2r  (r  -  degree  of  associa¬ 
tion  of  A)  and  r~=-l.  r 

s 

Using  expression  (6),  we  obtain 


'« -if = Dr  (t  -  0  !•  0  2  -dirf  -  y.A> (*•.  - *?> + 


(7) 


We  will  introduce  specific  designations  for  the  classes  of 
contact  sectors:  a  -  for  actyl  radicals,  p  -  for  phosphate  groups, 
m  -  for  phosphate  groups  in  which  at  least  part  of  the  ions  of 
hydrogen  are  replaced  by  the  ions  of  metal,  and  s  -  for  molecules 
of  carbon  tetrachloride.  Then  expressions  (.7)  acquires  the  form 

lg  »  f  (*  -  3)  2MT  I V'«  (U*"‘  “  V'P>  +  ~  U»]  + 

+  Vm(U . -  Ump)  -I-  [V,  -  1>  ( UM  -  U„) |. 
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it  follows  from  formula  (8)  that  the  ratio  of  activity  coef¬ 
ficients  being  considered  in  this  case  does  not  depend  on  r,  but 
evidently  depends  on  i  -  the  number  of  particles  M,  entering  into 
the  composition  of  the  given  complex,  which  should  be  taken  into 
account  during  summation. 

For  evaluating  the  interval  of  possible  changes  we 

1r 

will  use  the  following  simplifications:  U  =U-  — U  =co  and 

ms  ps  am  ap 

U  — -U-^=U  IL^.  The  degree  of  polymerization  of  mono-N-octyl- 

mm  pm  pm  pp  J 

phosphoric  acid  is  great  [24].  Therefore,  instead  of  (8)  we  get 


TMir  . _'  a  (s  -  3)  <» 

Ig — -».<»- -1 •  ssrftnt  (*-*). 
IT 


(9) 


where  -  total  concentration  of  A  in  the  organic  phase  and 


^  MM  flirt 

-  V,,  ■ 


Evaluations  of  the  variables,  analogous  to  C6  »  which  were 
obtained  in  a  number  of  works  [25-30],  showed  that  the  value  z* 
in  modulus  does  not  exceed  300  cal  •mole’"1.  On  the  other  hand,  in 
the  work  of  Tret'yakov  and  Rudakov  [31]  It  was  found  that  in  the 
most  diverse  liquids,  in  which  strong  hydrogen  bonds  are  absent, 
the  free  energy  of  intermolecular  interactions  lies  in  the  interval 
from  3  to  5  kcal •mole’*'5' .  This  interval  answers  to  a  change  in 
from  0.6  to  1.7. 

The  results  of  calculations  based  on  equation  (9)  with 

CO —~t  j  0 and  t=25°C  are  given  in  Table  1. 

It  follows  from  the  results  of  the  calculations  given  in  Table 
1  that  in  an  interpretation  of  experimental  data  based  on  the  dis¬ 
tribution  of  metals  in  extraction  systems  of  the  type  being  consid¬ 
ered  in  a  general  case  one  cannot  disregard  the  change  in  the  ratio 
YMir/y.r  with  a  concentration  of  extracting  reagent  C^>1  mole,l~“' 
and  with  i > 1 . 
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Table  1 

Greatest  possible  values  of  the  modulus  of  variable 


different  coordination  numbers  of  a  quasi-lattice  liquid 


for 


cn 

(K  M) 

Mo^yjik  If  JSE  npM  i  =  t 

_  _ Tr  _ _ 

t  =S*» 

*-=« 

taCD 

0.0! 

0.0001 

0.0003 

0.0006 

0.05 

0.0007 

0.0016  ‘ 

0.0028 

0.1 

0.0014 

0.0032 

0.0056 

0.5 

0.0070 

0.0158 

0.0280 

1.0 

0.0150 

0.0324 

0.0560 

Key:  Cl)  C.  (in  M);  (2)  Modulus  of  lg  when  1=1. 

A  if 


The  system  nitrate  of  tri-N^octylamine  (A)  -  solvate  (N)  - 
benzene  (S).  During  extraction  of  nitric  acid  by  an  amine  salt 
a  molecule  of  the  nitric  acid,  by  means  of  a  hydrogen  bona,  is 
connected  to  a  nitrate  group  of  the  amine  salt,  forming  a  solvate. 
Both  the  amine  salt  and  the  solvates  apparently  polymerized,  al¬ 
though  not  to  a  strong  degree  [1],  A  special  feature  of  this  sys¬ 
tem  is  that  the  dimensions  of  the  molecules  of  nitric  acid  and 
benzene  are  close  to  each  other.  Consequently,  the  Solvation  of 
amine  salt  by  nitric  acid  should  be  accompanied  by  growth  of  the 
Statistical  measurements  of  the  particle.  If  r,  just  as  previously, 
denotes  the  degree  of  polymerization  of  the  amine  salt  or  the  sol¬ 
vate,  then  r^r®5*r,  r,|ir=5*r+i  and  rg=l.  The  indices  a,  n  and  s 
denote  the  classes  of  the  contact  sectors  belonging  to  the  alkyl 
radical,  the  nitrate  group  and  the  molecule  of  benzene  respectively. 
In  this  case  equation  (6)  acquires  the  form 
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,(••;+  vm\\  -  2i' „  -  v,)  -  _  v,)\ } , 


(10) 


where  •»„  =  Vtj  - 1  (Vtt  +  v„). 
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Table  2 

Greatest  possible  values  of  the  modulus  of  for  different 

concentrations  of  amine  salt  (G^)  and  nitric  acid  (C^),  t=25°C 


(D 
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Mu/iyjit.  Ig  b 
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lif  JlpH  f  ss  I 
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I*  JW> 

cM/Ca 

ict 

l  . 

(l.l 

|  0.5 

|  1.0 

0.1 

1  0.5 

1  o’ 

_ _ 

2  =  00 

2  -+  8 

0.01 
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0.001 

0.001 

0.11*1 

0.05 

0.<X)3 

0.004 

0.004 

o.oou 

0.0015 

o.u* 

0.1 

0.007 

0.1X18 

0.009 

0.012 

0.014 

o.ur- 

0.5 

0.034 

0.038 

0,043 

o.ooo 

0,008 

o.o7: 

1.0 

0.007 

0.075 

0.08U 

0.119 

0.135 

o.lS. 

Key:  (1)  C,  (in  M);  (2)  Modulus  of 
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when  i=l. 


Keeping  in  mind  that  the  numerical  evaluations  are  carried 
out  for  a  concentration  of  reagent  up  to  i  mole ‘l-1,  the  following 
simplifications  should  be  recognized  as  permissible: 

ViVj^Vi,  and  .  Prom  equation  (10)  we  obtain 
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The  results  of  the  calculations  using  equation  (11)  are  given 
in  Table  2. 

As  it  follows  from  the  results  of  the  calculations  (Table  2), 
the  need  for  taking  into  account  the  change  in  the  coefficients  of 
activity  arises  here  at  lower  concentrations  of  reagent  than  in  the 
case  of  extraction  of  metals  by  mono-N-octylphosphdric  acid.  Con¬ 
sequently,  any  model,  intended  for  explaining  data  on  the  distri¬ 
bution  of  substances  in  systems  with  amine  salts  from  the  point  of 
view  of  chemical  events  in  the  organic  phase,  cannot  ignore  this 
fact . 
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Conclusions 

1.  A  variant  of  a  molecular-statistical  model  of  an  organic 
solution  (extraction  system)  is  proposed  which  includes  associated 
reagents  of  the  type  of  salts  of  alkylammonia  and  or  eranopho sphor ic 
acids.  The  statistico-thermodynamic  equations  take  into  account 
the  contribution  to  the  chemical  potential  of  a  particle  in  the 
mixture  of  the  specific  interactions,  dimensions  and  symmetry  of 
the  particles,  and  also  the  change  in  the  molecular  force  field 

of  the  medium  with  a  change  in  the  composition  of  the  system. 

2.  On  the  basis  of  equations  obtained  for  the  systems 
mono-N-octylphosphoric  acid— -complex — carbon  tetrachloride  and 
trioctylammonia  salt — nitric  acid — benzene  numerical  evaluations 
are  given  for  the  boundary  concentrations  where  it  is  necessary  to 
take  into  account  the  change  in  the  coefficients  of  activity  of 
the  separate  molecular  individuals  in  the  interpretation  of  data 


oh  distribution. 
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